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The properties of VLF propagating modes are briefly reviewed and simplified equations 
are presented which can be employed in calculating the fields produced. Experimentally deter- 
mined excitation factors are compared with theoretical curves obtained by Wait and found to 
agree rather closely. Attenuation rates are shown as a function of frequency and are found 
to agree rather closely with calculated values using a proposed simplified perturbation solu- 
tion for attenuation rates based on reflection coefficients at the ground and ionosphere. 
When employing values of ionospheric reflection coefficients recently presented by Wait and 
Walters, the nonreciprocal effects in attenuation rate are found both theoretically and 
experimentally to be much greater in the 10 kc/s region than in the 20 to 30 kc/s region. 
Finally, experimental field strength versus distance data are compared with mode calcula- 
tions and found to compare well all the way from 1 megameter out to and including fields at 
the antipode (20 mcgameters). 



1. Introduction 

The majority of VLF transmitting installations 
employ frequencies lying within the 10 to 30 kc/s 
region. Within this frequency region, and particu- 
larly at great distances, most of the energy can be 
considered as being propagated as a guided wave in 
the space between the earth and a conducting iono- 
sphere. Observed electromagnetic fields within the 
earth-ionosphere cavity indicate that the actual 
physical case can be well approximated with a model 
such as shown in figure 1 where the earth is considered 
as highly conducting with a reflection coefficient of 
+ 1 and the ionosphere is highly reflective with a 
reflection coefficient of approximately —1. For 
these idealized assumptions, the first order transverse 
magnetic (TMJ wave will have electrical field lines 
similar to those shown on the left hand side of the 
figure. The manner in which the magnitude of the 
vertical electric field, E z , and the horizontal electric 
field, E y , vary with height, z, is shown on the lower 
left hand side of the figure. The scale chosen for this 
figure is representative of what is expected with a fre- 
quency of 15 kc/s. On the right hand side of the figure 
are shown the electric field lines for the TM 2 ; i.e., 
the second order transverse magnetic wave. The 
amplitude-versus-height variations of the vertical 
electric and horizontal electric fields are again in- 
dicated at the bottom of the figure. The electric 
field lines are only shown within the cavity and no 
attempt is made to indicate the field lines lying 
outside the boundaries in the actual physical case 
where conductivities are not infinite. 

For distances where near field effects can be 
neglected, Wait [1962, p. 199] has shown that the 
vertical electric field E 2 can be represented in integral 
form or as the summation of an infinite number of 



modes. In the practical case for distances less than 
about 1 Mm, the total electric field at VLF can be 
well approximated as the sum of the groundwave and 
one or two sky reflected rays. This field is 



m=2 or 3 

E 2 =E l , g + £ E^ 

w = l 



(1.1) 



X/2<V/<1 Mm 



where E z<g is the groundwave contribution and E z%m 
is the mth hop sky wave contribution. 

For distances greater than 1 Mm, the vertical elec- 
tric field can now be represented as the sum of a 
relatively small number of waveguide modes l 



72 = 1 tO 3 

E^ S E,, n 

n = l 



(1.2) 



1 Mm<d. 



The field produced for a given mode can be written 
in decible terms as 

E z , n [db, v/m] = 10 Log P r +K' n - 10 Log/- 10 Log h it t 

+ 10 Log A Ut t — 10 Log (a sin d/a) 

-10 Log h itr +10 Log A n ,r+20 Log 

G h , n -L d , n -ad/W (1.3) 

where: E z , n [db, v/m] is the vertical electric field 
produced by the nth mode in db relative to 1 v/m, 



i At VLF the n=0 mode which is actually a TEM mode is highly attenuated 
and is not expected to contribute to the observed fields. 
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Figure 1. Electrical field structure for wave guide modes. 



P r is the power radiated into the half space above 
the earth in watts, K' n is a constant which relates 
field strengths at the surface to power radiated and 
is equal to 104.3 for modes of n=l,2- -, j is the 
frequency in cycles per second, h ut and h itT are the 
effective ionospheric heights at the transmitter and 
receiver locations, nominally 7X10 4 m for day and 
9 X 10 4 m at night. A n< t is the mode excitation modi- 
fication factor which is defined as the ratio of power 
launched into the concentric spherical shell guide rel- 
ative to that for a flat guide with perfectly con- 
ducting earth, a is the earths radius ^6.4 X10 6 m, d 
is the path distance, (a sin d/a) accounts for energy 
spreading in the spherical guide structure, A wr is the 
mode field modification factor for a curved earth at the 
ionospheric heights applicable at the receiver, 20 Log 
G^, n is the height gain factor in decibles relative to the 
field at the surface, i.e., z=0, L din is a factor in deci- 
bels which accounts for losses at discontinuities 
along the path, such as abrupt changes in surface 
conductivity, or in guide height such as can occur 
for a sunrise-sunset boundary, and a is the effective 
attenuation rate in db/Mm which is normally 
different for each mode number. 

As we will see later, this idealized model and 



associated equations is able to account rather closely 
for the observed variations in field strength for 
transmitters and receivers situated within the earth 
ionsophere cavity. 2 The important problems faced 
in obtaining numerical values for the factors in 
(1.3) consist primarily in determining: (1) the 
efficiency of a given transmitting installation in 
launching radiofrequency energy into the half space 
above the earth, (2) the efficiency with which this 
energy couples into the propagating mode, (3) the 
manner in which the energy spreads throughout the 
concentric spherical shell cavity, (4) the height gain 
characteristic, (5) the discontinuity loss factor, and 
(6) the attenuation rate in the last term which 
accounts for losses along the path. 

The effective radiated power P r can be determined 
from near field measurements applied to well-known 
relations. The frequency, of course, is known and 
the other geometrical factors involving distance, 
height to the ionosphere, and the spreading terms 
can be rather precisely determined. The principal re- 



2 Normally in the 1 to 10 Mm distance region the fields are considered as being 
produced primarily by the TM2 and TMi modes. The possibility of TE modes 
contributing must not be neglected particularly at short or very long distances 



maining unknowns relate to the excitation efficiency 
and at t (Mi t nation rates. 

In order to determine attenuation rates, Wait 
[1958] employed a technique in which he plotted 
the sum of the vertical electrical field in decibels 
relative to 1 v per meter and the spreading term 
10 Log {a sin d/a), as a function of the distance. 
This technique has been employed by the authors 
for a large number of observed field strength versus 
distance contours to obtain both excitation factors 
and attenuation rates. A typical E z (a sin d/a) 
versus d curve is shown in figure 2 for the trans- 
mitting station of LCM of Stavanger, Norway 
operating at a frequency of 24.4 kc/s. The data 
was taken by Round et al. [1925]. For a distance of 
1 km, the energy spreading term (a sin d/a) becomes 
equal to zero and if the radiated power, frequency, 
and ionospheric height are known, it is possible 
to arrive at a 1 km intercept point which for all 
practical purposes can be plotted for d=0 Mm 
(megameters) on the curve in figure 2. For this 
particular station, the calculated (/={) intercept 
value is 89.5 decibels while the actual intercept is 
seen to be some 9 db below the flat earth calculation. 
Since the attenuation is negligible for the 1 km 
distance, the resulting field strength difference is 
actually equal to 20 Log A which is the total modifi- 
cation required for the spherical earth case; i.e., 
20 Log A = 10 Log A f +10 Log A r if the ionosphere 
height is the same at transmitter and receiver. 
The slope of the curve, since the distance spreading 
term has been added to the observed fields, is seen 
to be constant and for this particular case is equal 
to about 2.3 db/Min. The attenuation rate, a, 
is determined by losses at the earth and ionosphere 
boundaries and will be shown later to be dependent 
upon many factors including: frequency, diurnal 
effects, seasonal effects, latitude effects, solar activ- 



ity, meteoric activity, direction of propagation 
relative to the earth's magnetic field, and earth's 
surface conductivity. 

2. Excitation Factors and Height Gain 

The power radiated into the half space above the 
ground surrounding the antenna is given by the 
well-known relation 



Pr, 



half space" 



--P 1607r 2 (A e /X) 2 



(2.1) 



where / is the rms antenna current, h e is the effec- 
tive antenna height in meters, and A is the wave- 
length in meters. When the antenna is considered 
as being located between two parallel plates, the 
fields at a large distance from the source can be 
considered as though they were produced by a set 
of images producing an equivalent line current as 
shown by Wait [I960]. Employing the field pro- 
duced by this equivalent line source, it can be shown 
that- the power launched into a given spherical earth 
mode relative to the total half space radiated power is 



P 

■*- r, n = 



(0,1.2.. .)_ 



77 () XA 1.12X10 8 A 



Pr, 



half space 



:',2(Wr 



fh 



(2.2) 



where rj is the impedance in free space, / is the fre- 
quency in cycles per second, h is the ionospheric 
height in meters, and A is the mode excitation mod- 
ification factor. For a spherical earth, the variation 
in electromagnetic field strength with height changes 
from that of the flat earth case. 'Phis change in 
field strength with height modifies, by a factor A, 
the amount of energy which is coupled into the 
spherical cavity or concentric spherical cavity. In 
addition, the field strength observed at the surface 
by the receiver is reduced by a factor A 1 ^ for a fixed 
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Figure 2. Typical E (a sin d/a) versus d curve. 



energy flow in the guide. The result, assuming h t 
is constant, is that the vertical electric field calcu- 
lated for a flat earth must be modified by A to the 
first power. Wait [1962] has derived values for the 
modification factor A which are shown in figure 3 
for a typical day time ionospheric height of 70 km. 
The curves are for the first and second order modes 
where the earth and ionosphere boundaries are 
considered as perfectly reflecting. 

In order to determine A experimentally, a large 
number of field strength versus distance curves were 
plotted (in a manner similar to that shown in fig. 2) 
and the resulting experimental data points shown in 
figure 3 appear to group around the n=l curve over 
quite a large frequency region. The few data points 
shown at the upper right hand side of the figure 
near the n = 2 curve, were obtained from line inter- 
cepts which appeared to be representative of second 
order mode propagation. These line intercepts cor- 
respond to a rather high attenuation rate portion 
of the field strength versus distance curve that in 
general is only applicable out to several mega- 
meters from the source. As might be anticipated, 
these effects are seen essentially only at the higher 
frequencies; i.e., about 20 kc/s and up. 

It is interesting to observe that on several occa- 
sions the excitation factors appear to lie appreciably 
below the theoretical curve. For example, at 21 
kc/s on one occasion the excitation factor was near 
the curve while for the other example shown, the 
excitation modification apparently produced was 
about 5 db below the theoretical value. The par- 
ticular field strength versus distance curve associated 
with this low excitation factor had a rather low 
attenuation rate, and it is likely that the propaga- 
tion conditions were approaching the earth detached 
case as described by Wait [1962] where most of the 
energy is concentrated near the ionosphere and any 
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Figure 3. Excitation factors (theoretical and experimental). 



coupling to antennas at the earth's surface is very 
small. Obviously, in this case the vertical electric 
field versus height will not be a maximum at the 
earth's surface. 

In the experimental data shown, an attempt was 
made to group the examples into propagation in 
both the north-south direction as compared to that 
in the east and the west directions. In addition, 
those paths which were launched over sea water 
and others over land were also examined. It was 
impossible from the rather limited data available 
to determine for certain any trend such as an 
increased excitation loss A for launching over land 
or in any specific direction. Since, as we will see in 
the following section on attenuation rates that the 
ionospheric reflection coefficients are different for 
various directions of propagation, it is anticipated 
that a close examination of excitation factor may 
indicate some sensitivity in A to direction of propaga- 
tion relative to the earth's magnetic field, but at 
present it must be concluded that we have not 
demonstrated this point. It is interesting to observe 
here that Wait and Spies [1963] in their recent 
descriptions of height gain for finitely conducting 
surfaces, show that the excitation factor shows a 
somewhat lower reduction, i.e., more efficient excita- 
tion, for poorly conducting earth than for perfectly 
conducting earth. This of course is for a condition 
where the whole earth's surface is uniformly poorly 
conducting. For the case where the guide parameters 
are primarily determined by a sea-ionosphere set 
of boundaries and the transmitter is located some- 
what inland, there may be a cutback modification 
required which is expected to show a decrease in 
excitation efficiency when the antenna is situated 
over a finite amount of poorly conducting terrain. 
For normal conducting grounds at VLF this modifica- 
tion is expected to be rather slight; however, as 
frequency increases, particularly up into the LF 
band, this type of modification can be very pro- 
nounced as shown for arctic paths by Watt, Maxwell, 
and Whelan [1959]. 



3. Attenuation Rates 

The effective attenuation rate used in the wave- 
guide modal equation represents an average loss 
per unit distance for a discrete propagation path 
and time; therefore, any change in the electrical 
properties of the ionosphere with time and any 
change in the electrical properties of the earth 
and ionosphere over various paths will give a cor- 
responding change in the effective attenuation rate. 
It is well known from work by Budden [1952], 
Barber and Crombie [1959], and Wait [1961] that 
the attenuation rate is also affected by the earth's 
magnetic field. This effect, which is produced by 
the modified electron motion in the presence of a 
magnetic field, may be treated as an earth's magnetic 
field factor M(<j> a ) 3 which modifies the attenuation 



3 Wait uses the symbol P for this factor. 



rate as follows 

■jir/, n _ Attenuation with Magnetic Field . , 
^^"Attenuation without Magnetic Field ^ } 

whore 

</> a = azimuthal direction of propagation with 
respect to the horizontal component of 
the earth's magnetic field. 

The value of M is also dependent upon the 
geomagnetic latitude which determines the magnitude 
of the horizontal components of the earth's field. 
Since M is the influence of the earth's magnetic 
field on VLF propagation by virtue of modification 
of the reflection coefficients at the ionosphere, any 
additional loss in the waveguide produced at the 
lower, i.e., earth's surface, boundary is not dependent 
upon M. 

To a first approximation, the effective loss over 
a given propagation path d f megameters long with 
various effective ionosphere and ground conductivities 
along the path can be considered as the sum of 
all the incremental losses along the guide surfaces. 
The effective loss can be written as 

Ni """ N e 

ad/lO«=ad' =52M k (<t> a )a tt *A4 +2 ^ g , ,A^ (3.2) 

/;=1 ; = 1 

where: N t and N g are the number of ionospheric 
and ground path increments considered, M k 
{<t>a)a ltk Ad k is the loss contributed by the ionosphere 
over the incremental distance Ad' k , and Aa g<j Adj is 
the loss contributed by the earth over the incre- 
mental distance Ad/. For paths where the boundary 
surfaces can be considered as uniform, the attenua- 
tion rate becomes 



a=M(<t> a )at+Aag t 



(3.3) 



The contribution due to ground loss can be shown 
to be related to ground conductivity by the ex- 
pression obtainable from Wait [1957] 



Aa,[db/Mm]--46/ /2 o-- 1/2 A- 1 ; 



(3.4) 



where : / is the frequency in cycles/second, a- is the 
ground conductivity in mhos/meter, and h is the 
effective ionospheric height in meters. This ex- 
pression is valid only when / is low enough that the 
mode is not becoming earth detached; i.e., about 
15 kc/s or less for n=l. The conductivity value 
must also not be too low for this relation to apply. 
The earth detached mode may be seen by the more 
rigorous computations of Aa g by Spies and Wait 
[1961] in which a maximum value of Aa g is reached 
with an increase of frequency, and a further increase 
will decrease Aa g . 

A great deal of work has been conducted in cal- 
culating the attenuation rates of the ionosphere 
and earth of different conductivities by a number of 
authors. The results are presented either in the 



form of attenuation rate or in the form of reflection 
coefficients at the earth and ionosphere. 

The attenuation rate can be obtained from the 
reflection coefficients as a perturbation solution which 
assumes that 



a[db/Mm]: 



g,[db]+gjdb] 
2d[ 



(3.5) 



where: R g and R t are the reflection coefficients of 
the ground and ionosphere respectively, and 2d[ 
is the average distance in megameters containing one 
ground and one ionosphere reflection at mode res- 
onance. An approximate solution for d r that is 
valid over an appreciable range of frequencies yields 

a[db/Mm]~1.2X 10 13 (JS,[db] 

+R t [db]) (2im-<l> g -fa)/fli 2 (3.6) 

where 4> g and <t> t are the phase lag upon reflection at 
the ground and ionosphere, respectively. 

Values of R { from Wait and Walters [1963] are 
shown in figure 4 as a function of frequency and 
parametric in the term I which is the conductivity 
gradient of the ionosphere. 4 Typical conductivity 
gradients are shown in figure 5 for both day and night. 
The conductivity gradient term, Z, may be calculated 
from these curves by the following relationship: 



h 



'2.3Log(<r/tr ) 



(3.7) 



where : and : are the ionosphere heights correspond- 
ing to the conductivities <r Q and a, respectively. 

The ionospheric conductivity neglecting the earth's 
magnetic field is given by Jordan [1950] as 



Ne 2 



m(i> 2 +co 2 ) 



~2.78X10- 8 iV> 



(3.8) 



where: e and m are the electron charge and mass, 
N is the number of electrons per cubic meter, and 
v is the effective collisional frequency in collisions 
per second. Typical daytime and nighttime values 
of the conductivity gradient are Z day ~2.6 km and 
4iight^l.5 km. 

Anticipated sea water path attenuation rates can 
be calculated with the aid of (3.6) assuming that 
/=20 kc/s, R g ^0 db, <t> g =0, h a&y =70 km, and 
^ni g ht=90 km. From figure 4, i? iday ~4.5 db and 
#*night=3.5 db, and from Wait and Walters [1963], 
<t>i day =2.4 and fa night=2.7. The resulting calculated 
attenuation rates are 

a day =2.2 db/Mm 
^night^O.9 db/Mm. 



4 Note that Wait uses the parameter £ which equals 1/Z. One of the most 
significant recent contrihutions to VLF propagation theory is contained in the 
paper, "Reflection of electromagnetic waves from a lossy magneto plasma," 
by J. R. Wait and L. C. Walters [1963a]. The model employed in this paper is 
a close representation of what we presently expect for the ionosphere. The re- 
sultant values of R% now reflect the influence of the earth's magnetic field as well 
as finite conductivity gradients. 
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Figure 4. Reflection coefficients versus frequency for nonsharply bounded ionosphere. 
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Figure 5. Typical ionosphere conductivity gradients. 



Curves of experimentally determined daytime 
values of attenuation rates are shown in figure 6 
as a function of frequency and parametric in direc- 
tion of propagation. Much of the data for this 
figure came from curves similar to figure 2. The 
curves were drawn through the average of all the 
points within a plus or minus 45° sector of the given 
direction. 

For frequencies above 25 kc/s, existing data is 
contaminated by higher order mode contributions. 
Because of this, the dashed portions of the curves 



are based primarily on theoretical values by Wait 
and Walters [1963a]. 

Wait [1964] has calculated the change in attenua- 
tion rate produced by an exponential ionosphere 
relative to that previously calculated using a sharply 
bounded ionosphere [Spies and Wait, 1961]. The 
resulting exact solutions of a for an exponential 
ionosphere (h=70 km day and 90 km night) may 
be seen in figures 7 and 8 for the first and second 
order mode respectively. The curves were obtained 
by using typical values of the conductivity gradient 
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Figure 6. Attenuation rates for predominately sea ivater 
paths versus frequency for moderate and low sunspot numbers 
(experimental) . 
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Figure 7. Attenuation rates over predominately sea water 
paths for n = l (experimental and theoretical) . 

expected for day and night. Indicated on the 
curves are experimentally determined values of 
attenuation rate. It can be seen that the theoretical 
values are in good agreement with those obtained by 
experimentally collected data for the first order 
mode. For the second order mode, the experi- 
mental values obtained are much more uncertain 
which may account for the only fair agreement 
during the night. This uncertainty may be due to 
the randomness of nighttime propagation condi- 
tions; i.e., variable I and h. In addition, the distance 
involved in which the second order mode is pre- 
dominate is usually very short for accurate a deter- 
minations. 

4. Field Strength Versus Distance 

Figure 9 shows the observed vertical electrical 
field strength as a function of frequency for a typical 
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Figure 8. Attenuation rates over predominately sea water 
paths for n = 2 (experimental and theoretical) . 



west to east propagation path. Also shown or the 
figure are the calculated ground wave and modal 
fields using values of the excitation factor A as shown 
in figure 3 and expected values of attenuation as 
shown for n=l in figure 6 and n=2 in figure 8. At 
distances less than 300 km, the groundwave is pre- 
dominate and is practically independent of frequency 
over the VLF range. This may be seen from the 
matching of data points at 20 and 16 kc/s. At ap- 
proximately 100 km, an interference pattern begins 
to appear which is likely due to interference between 
the groundwave and the first hop sky wave. Near 
600 km, the groundwave and sky wave are approx- 
imately of equal amplitude. Since their phase versus 
distance rates are different, a strong cancellation, or 
null, is expected in this region. 

Beyond 1 Mm where we are now considering mode 
theory, the first order mode is apparently predomi- 
nate and has an attenuation rate of approximately 
2 db/Mm. Beyond about 16 Mm, it can be seen that 
the converging energy due to antipode focusing actu- 
ally increases the fields at a greater rate than the loss 
due to attenuation. This buildup in field strength 
continues until the antipode is reached. It can be 
seen that the observed antipode value [Garner, 1962] 
is very close to the calculated value. Included on 
the figure are the unattenuated mode fields for the 
first and the second order modes. 

If the transmitted frequency is increased, the 
attenuation of the second order mode will decrease 
rapidly and at the same time the excitation of the 
first order mode will decrease. It then is apparent 
that as the frequency increases the second order mode 
will predominate to greater ranges. Figure 10 is a 
curve of vertical field strength as a function of fre- 
quency for east to west propagation for data in the 
21 to 25 kc/s region. The excitation factors were ob- 
tained from figure 3 and the attenuation rates were 
chosen to best fit the data. The first large null in 
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field strength is now seen to occur at 700 km and is 
again likely caused by cancellation between the 
ground wave and the first hop sky wave. 

Considering mode theory for distances greater 
than 1 Mm, the second order mode appears to be the 
predominate mode out to about 3 Mm where inter- 
ference now occurs between the first and second 
order modes. Beyond 4 Mm the first order mode 
appears to be predominate. On some occasions, the 
second order mode has been observed to be predomi- 
nate at nighttime for distances up to 6 Mm or more 
for frequencies above 18 kc/s. 
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